Logical and Arithmetic Circuits in Belousov Zhabotinsky Encapsulated Discs 
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Excitation waves on a sub-excitable Belousov Zhabotinsky (BZ) substrate can be manipulated 
by chemical variations in the substrate and by interactions with other waves. Symbolic assignment 
and interpretation of wave dynamics can be used to perform logical and arithmetic computations. 
We present chemical analogs of elementary logic and arithmetic circuits created entirely from inter- 
connected arrangements of individual BZ encapsulated cell like discs. Inter-disc wave migration is 
confined in carefully positioned connecting pores. This connection limits wave expansion and unifies 
the input-output characteristic of the discs. Circuit designs derived from numeric simulations are 
optically encoded onto a homogeneous photo-sensitive BZ substrate. 

PACS numbers: 82.40.Ck, 89.20.Ff, 89.75.Fb, 89.75.Kd 



I. INTRODUCTION 

Unconventional Computing [1] is a field of study ded- 
icated to exploring alternate computational strategies, 
structures and mediums contrasting contemporary digi- 
tal von Neumann architecture machines [2]. Innate se- 
rial execution and storage processing of these machines 
presents profound throughput and energy restrictions, 
limiting practical application of some classes of algo- 
rithms [3]. Conversely computational capable organic 
adaptive systems are innately parallel and distributed. 
Highly parallel and distributed computation can be sup- 
ported in chemical reaction-diffusion (RD) systems (J] 
such as the Belousov Zhabotinsky (BZ) reaction [5]. One 
bench mark of these systems is the ability to replicate 
components used in conventional computation, such as 
logic and arithmetic circuits. Logic gates are the phys- 
ical embodiment of Boolean logic operations that form 
the foundation for digital computation. Circuits of logic 
gates can be connected to create machines capable of 
performing ^Universal Computation' [6[. 

The chemical implementation of logic gates were first 
described by Rossler sustained in a bistable chemical 
medium [7] . Logic gates have also been demonstrated us- 
ing chemical kinetics [8-11]. The binary value assignment 
to the presence or absence of individual RD waves modu- 
lated by the substrate geometry and other waves presents 
an alternative approach. Toth first described logic gates 
with BZ waves along these lines [3, [l3|. In that system 
RD waves travelled along thin capillary tubes. These 
basic principles have since led to numerous adaptations, 
for example Steinbock implemented logic gates by 'print- 
ing' a catalyst of the BZ reaction onto a facilitating 
medium In a simulated study Motoike presented 

(amongst other designs) logic gates reliant upon geomet- 
ric patterns of passive diffusion boundaries embedded on 
an excitable field ^5|. Sielewiesiuk also created logical 
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functions from using passive diffusion lines where func- 
tion arose at the intersection of perpendicular wave chan- 
nels [16]. Adamatzky presented alternative logic gate de- 
signs by combining the principles of collision based com- 
puting on an precipitating chemical substrate and 
in- vitro 1 191 . Similar or derivations on the above studies 
such as: [20| using the FitzHugh-Nagumo model [^T, [22| 
reaction, or [23|] co-evolutionary system, or cardiac cell 
model on a GPU [IJ and amongst others [25] have all 
implemented logic functions or gates either in simulation, 
in-vitro or both. 

Current experiments with lipid coated vesicles of ex- 
citable media at our (and our partners [26^) laboratories 
have led to a novel geometric BZ mediated approach to 
creating logic and arithmetic functions. Representative 
of a 2 dimensional approximation of BZ vesicles, networks 
of interconnected discs containing a sub-excitable BZ for- 
mulation can be arranged to create various circuits. Logic 
& arithmetic circuits and polymorphic gates have been 
simulated in homogeneous hexagonal and orthogonal net- 
works [2?, ^28^ and non- homogeneous arrangements [30| . 

In these simulations, geometric arrangements of discs 
are joined together by small connecting pores. The re- 
sulting pattern modulates a homogeneous BZ Oregona- 
tor simulation grid that has adapted for monochromatic 
photo-sensitivity [s^. The disc boundaries are created 
from rings of high intensity light to emulate an impen- 
etrable exterior membrane. Lower intensity light forms 
the interior disc area and connecting pores adjusted to 
support travelling waves. 

Logically symbolic waves are able to traverse the net- 
work modulated by interaction with pathways and other 
waves. The disc interior can be exploited for free space 
collision style reactions [17] whereas the pore loci and ef- 
ficiency can compartmentalise the resulting reaction [2^ • 
Circuits have been created from logical sub assembles in 
orthogonal and hexagonal networks [l^, [sOj. Function 
density can be increased by circuits that include varia- 
tions in relative disc size, pore efficiency and connection 
angles ^] . Using this later technique we have assembled 
a compact elementary arithmetic circuit. 
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Conceivably other cell geometries (such a grids of 
square units) could compartmentalise the cell function 
similar to discs. Nevertheless, the curvature of the wave 
front in combination with the disc geometry provides 
a natural and convenient unit and network structure. 
In comparison to our previous free space [32| or chan- 
nelled logic optically projected circuits [SS^ the combina- 
tion of the natural positively expanding wave curvature 
the disc geometry promotes functional inter-disc stability 
and connectivity. 

Use of a photo-sensitive adapted version of the Orego- 
nator model permits a simple migration from simulation 
to experiment. Circuit designs from the simulation can 
be projected directly onto an actual photo-sensitive BZ 
medium. Furthermore, the same substrate is capable of 
supporting different successive circuits. Detected activity 
in assigned output discs can trigger the projection of new 
alternate circuits onto the same substrate area. Those 
triggering waves can be captured into the new circuit, ef- 
fectively becoming inputs into the new circuit. Migration 
of our stabilising disc centric circuit simulations onto a 
reusable chemical substrate is the primary subject of this 
report. 

II. METHODS 

Designs are initially created in simulation and then the 
resulting geometry is then projected onto the BZ sub- 
strate. This substrate is a photo sensitive BZ system [34] 
with a Ru(bpy)|^ catalyst immobilised on silica gel j35j . 
The medium is oscillatory in the dark with composition 
[36]. Global illumination levels are manually adjusted to 
promote marginal excitation waves. Simultaneous mul- 
tiple wave initiations are optically instigated. Removing 
all inhibiting light promotes mass wave excitation across 
the BZ substrate. Normal light levels and experimental 
pattern are then restored with the exception of small cir- 
cular zones at desired initiation points. Wave fragments 
continue to emerge from these dark zones out into the 
marginal zone of the disc body. 

A. Numerical simulation 

Numerical simulations are based on a 2 variable ver- 
sion of the Oregonator model [37] as a model of the BZ 
reaction [5|, [38| adapted for photo-sensitive modulation 
of the Ru-catalysed reaction j3ll |. 




Variables u and v are the local instantaneous dimen- 
sionless concentrations of the bromous acid autocatalyst 
activator HBr02 and the oxidised form of the catalyst 



inhibitor Ru(bpy)3^. ^ symbolises the rate of bromide 
production proportional to applied light intensity. Bro- 
mide Br is an inhibitor of the Ru-catalysed reaction, 
therefore excitation can be modulated by light intensity; 
high intensity light inhibits the reaction. Dependant on 
the rate constant and reagent concentration e represents 
the ratio of the time scales of the two variables u and 
V. q is di scaling factor dependent on the reaction rates 
alone. The diffusion coefficients Du and Dy of u and v 
were set to unity and zero respectively. The coefficient 
Dy is set to zero because it is assumed that the diffusion 
of the catalyst is limited. Inputs where instigated by per- 
turbing the activator = 1.0 centrally in small circular 
disc areas, rad = 2 simulation points (SP). 

Parameter Value Description 



€ 


0.022 


Ratio of time scale for variables u and v 




0.0002 Propagation scaling factor 


f 


1.4 


Stoichiometric coefficient 






Excitability level (proportional to light level) 


u 




Activator HBr02 


V 




Inhibitor Ru(bpy)3^+ 


Du 


1.0 


Activator diffusion coefficient 


Dy 





Inhibitor diffusion coefficient 


Ax 


0.25 


Spatial step 


At 


0.001 


Time step 



TABLE I. Kinetic and numerical values used in numerical 
simulations of (Eq. [T|). ★ Varies between two levels, sub- 
excited (LI) and inhibited (L2), 0li = 0.076, 0l2 = 0.209 

Numerical simulations were achieved by integrating 
the equations using the Euler-alternating direction im- 
plicit (ADI) method, method [39] with a time step 
St = 0.001 and a spatial step 6 = 0.25. Experimental 
parameters are given in Tab. HI 

Networks of discs where created by mapping 2 dif- 
ferent (j) values (proportional to light intensity) onto a 
rectangle of homogeneous simulation substrate. To im- 
prove simulation performance the rectangle size was au- 
tomatically adapted depending on the size of the net- 
work, but the simulation point (SP) density remained 
constant throughout. The excitation levels, L1^L2 re- 
late to the partially active disc interiors and non-active 
substrate. 

Discs are always separated by a single simulation point 
(1 SP) wide boundary layer. Connection pores between 
discs are created by superimposing another small link disc 
at the point of connection (typically a 2^6 SP radius), 
simulation points have a 1:1 mapping with on screen 
pixels. The reagent concentrations are represented by 
a red and blue colour mapping; the activator, u is pro- 
portional to red level and inhibitor, v proportional to 
blue. The colour graduation is automatically calibrated 
to minimum and maximum levels of concentration over 
the simulation matrix. The background illumination is 
mono-chromatically calibrated in the same fashion pro- 



portional to 0, white areas are inhibitory and dark areas 
excited (for monochromatic reproductions, blue (v) has 
been suppressed and red {u) appears as hght gray, colour 
versions online). 

Wave fragment flow is represented by a series of super- 
imposed time lapse images (unless stated otherwise), the 
time lapse is 50 simulation steps. To improve clarity, only 
the activator (u) wave front progression is recorded. Fig- 
ure [T] illustrates the same wave fragment in both colour 
map {u & v) and final composite time lapse version (u). 




(a) (b) (c) (d) (e) (f) 



FIG. 1. (Colour online) Example the time lapse image cre- 
ation. The image shown in (f) is the accumulation of suc- 
cessive images shown from central wave initiation in (a) to 
extinction in (e). Time lapses periods are 50 time steps and 
the refractory tail of the inhibitor (u) shown in blue (gray 
inner shadow of propagating ring [(a)^(e)]) is not shown in 
the time lapse image (f) to improve clarity. 

Inputs are created by perturbing a small circular area 
of the activator (u) set to a value of 1.0 with a radius of 
2 SP in the center of the disc. In addition to explicit la- 
belling, all discs representing inputs and outputs are also 
highlighted with a blue and green border respectively. 

In a previous study we have shown that logic circuits 
can be created with uniform discs arranged in hexagonal 
networks [2^ , hexagonal packing being the most efficient 
method of sphere (disc) packing. Further opportuni- 
ties to modulate wave fragment behaviour are presented 
when disc size, connection angle and connection efficacy 
are combined in non- homogeneous networks. Disc size 
can be adjusted to permit or restrict internal wave in- 
teractions, producing either larger reaction vessel discs 
or smaller communications discs (Fig. [2^). Connection 
angle between discs can be used to direct wave colli- 
sions (Fig. [JJd) and connection efficiency can effect the 
wave focus (Fig. [2t). 



B. Experimental 

A Sanyo proxtrax multiverse projector was used to 
project the design outline of several simulated circuit 
designs from [30]. Wave activity was captured using a 
Lumenera infinity 2 USB 2.0 scientific digital camera. 
The open reactor was surrounded by a water jacket ther- 
mostatted at 20° C. Peristaltic pumps were used to pump 
the reaction solution into the reactor and remove the ef- 
fiuent (Fig.[3l). 

Sodium bromate, sodium bromide, malonic acid, 
sulphuric acid, tris(bipyridyl) ruthenium(II) chloride, 
27% sodium silicate solution stabilised in 4.9M sodium 
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(a) Relative disc (b) Connection angle (c) Connection 
sizes efficiency 

FIG. 2. (Colour online) (a) Discs as reaction vesicles in the 
LHS network or communication channels in the RHS network, 
(b) The effect of connection angle, the two input signals com- 
bine in the left network to produce an output. Adjusting 
the angle of the lower input in the right network alters the 
result of the collision and no output is produced, (c) The 
effect of connection efficiency. Large aperture (6 SP) connec- 
tion in the left network results in a broad spreading beam. 
Conversely a smaller (4 SP) aperture connection in the right 
network creates a narrow beam wave. 



hydroxide [40] and used as received unless stated other- 
wise. Ru(bpy)|^ was recrystahsed from the choloride 
salt with sulphuric acid. To create the catalyst loaded 
gels a thin layer chromatography pre-coated plates silica 
gel with 254nm ffuorescent indicator on glass was cut into 
5cm X 5cm and placed in a solution of 0.9m/ of 0.025M 
Ru(bpy)|^ and 12m/ of deionised water in a Petri dish 
for 12 hours. Gels were washed in deionised water to 
remove by products before use. 

The catalyst free reaction mixture was freshly prepared 
in a 30m/ continuously- fed stirred tank reactor (CSTR), 
which involved the in situ synthesis of stoichiometric bro- 
momalonic acid from malonic acid and bromine gener- 
ated from the partial reduction of sodium bromate. This 
CSTR in turn continuously fed a thermostatted open 
reactor with fresh catalyst-free BZ solution in order to 
maintain a nonequlibrium state. The ffnal composition 
of the catalyst-free reaction solution in the reactor was: 
0.42M sodium bromate, 0.19M malonic acid, 0.64M sul- 
phuric acid and 0.1 IM bromide. The residence time was 
30 minutes. 

The spatially distributed excitable ffeld on the surface 
of the gel was achieved by the projection of the disc pat- 
tern from the simulations. The light intensity of discs 
was controlled by computer. The pattern was projected 
onto the catalyst loaded gel via a data projector. Every 
10 seconds, a light level of 5.7mW/cm^ for 10ms during 
which time an image of the BZ fragments on the gel was 
captured. The purpose of this was to allow activity on 
the gel to be more visible to the camera. Captured im- 
ages were processed to identify chemical activity. This 
was done by differencing successive images to create a 
black and white image. The images were cropped and 
layered to show progression of a single image ffnally the 
disc boundary was superimposed on the images to aid 
analysis of the results. 
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FIG. 3. Experimental setup: A light sensitive catalyst 
Ru(bpy)|^ loaded silica gel is immersed in catalyst free BZ 
reaction solution in a thermostatted (G) Petri dish (E). Peri- 
staltic pump (F). Reactor with thermostatted reaction solu- 
tion and to remove effluent (J). The reaction solution reservoir 
(H) is kept in an ice bath. The heterogeneous network on the 
surface of the gel (D) is constructed by the projection (B) 
of a disc pattern generated by a computer (A) via a mirror 
assembly (C). Images captured with a Digital camera fitted 
with a 455nm narrow bandpass interference filter (I). 



III. RESULTS 

We report on a selection of experiments that success- 
fully migrated from simulation; a diode, NAND gate, XOR 
gate and an simple arithmetic adder circuit [30]. In the 
following images the inverse monochromatic wave pro- 
gression (left column) is illustrated adjacent to the re- 
spective colour simulation analogue (right column). Al- 
though the parameterisation of the simulation creates dif- 
ferent wave dynamics the projected geometry in each case 
is proportionally identical resulting in the same function- 
ality. 



A. Diode 

The diode junction constructed with BZ discs oper- 
ates by simultaneously exploiting a right angle junction 
and asymmetric pore size (Fig. |4]). A wave fragment (y) 
travels horizontally from right ^left (Fig. 2^). Wave frag- 
ments cannot survive when the fragment size drops below 
some critical level [41]. A constricted pore of (4 SP) at 
the perpendicular junction reduces the wave size close 
to termination in the proceeding disc. Wave develop- 
ment is momentarily marginal before slowly recovering. 
Delay in wave development prohibits the fragment en- 
tering the vertical connection (Movies Mov. la&b [42]). 
Conversely the opposing vertical fragment {y) enters the 
junction through a broad pore (6 SP). The resultant 
increased junction wave seed permits the fragment to 
develop more rapidly and dissipate into the horizontal 
discs line (Fig. [Dd). Therefore signal propagation only 
progresses from vertical to horizontal {x^y). Discs in 
the simulation have a diameter of 56 SP which results 
in a projected gel disc of 10mm . This diode junction 
has been used to inject a unidirectional wave fragments 



into a type of simulated 1 bit memory where unidirec- 
tional waves indefinitely circulate around a loop of discs 
until reagent exhaustion or collision with an opposing 
wave [s^. 




FIG. 4. (Colour online) Diode junction {x^y = 0, y^x = 1). 
Left column sequences show wave (signal) development com- 
pared alongside analogous numerical simulations on the right, 
(a) Signal propagates horizontally (x^y). A narrow band 
pore (4 SP) at the 2'^^ (right angle) junction prohibits prop- 
agation downwards to y. (b) Signal propagates vertically 
(y^x). A broad band pore (6 SP) at the 2^^ (right an- 
gle) junction connection permits the signal to expand hor- 
izontally towards x. Simulation to gel projection ratio is 
56 SP : 10mm (Mov. la&b [42]). 



B. NAND logic gate 

A NAND gate design is shown in Fig. [5] and animated 
in movies Mov. 2a^c [42|. The gate is a conjunction 
of an AND gate and an inverter or not gate. The not 
gates operates horizontally in the bottom four discs. A 
logical TRUTH or '1' (supply wave) is always presented 
simultaneously with the inputs x h y. In the case of 
(x, y) = (0, 0) a wave propagates horizontally unimpeded 
terminating in the output disc {z = 1) (Fig. [5^). Pre- 
sentation of input cases (x, ?/) = (1, 0) | (0, 1) also have no 
effect on the progression of the supply wave and in both 
cases z = 1 (Fig.fSb). In the final instance (x^y) = (1, 1) 
the inputs collide in the central disc resulting in a per- 
pendicular ejection wave. The downward ejection wave 
fragment propagates into the horizontal not line below 
and deflects and causes the extinction of the supply wave 
leading to a negated output {z = 0) (Fig. [St). The NAND 
gate is significant because of its universal applicability 
in the construction of all other logic gates. This design 
is illustrative of the modular facilitation of BZ disc net- 
works creating a NAND gate from the conjunction of an 
AND and a not gate. 




FIG. 5. (Colour online) Two input NAND gate {z = x • y). 
Left column sequences show wave (signal) development com- 
pared alongside numerical simulations on the right. (a) 
z = l,{x,y) = [0,0]. (b) z = 0,ix,y) = [0,1]|[1,0]. (c) 
^ = 0, (x, y) = [1,1]- Pores of 6 SP interconnect 56 SP0 discs 
with simulation to gel projection ratio of 56 SP : 10mm 
(Mov. 2a^c HI). 



C. XOR logic gate 

Figure. [6] demonstrates the XOR function. As with 
the proceeding NAND gate, a logical truth ('!') is sup- 
plied in two discs (centre top and bottom right) in syn- 
chronisation with the two inputs {x & y). Presenta- 
tion of (x, y) = (0, 0) allows the vertical supply wave 
to deflect the horizontal supply wave resulting in logi- 
cal false output z = (Fig. Oa). Complimentary inputs 
(x^y) = (0,1) & (1,0) deflect the vertically travelling 
supply wave allowing the horizontal supply wave reach- 
ing the output disc {z = 1) (Fig. Ob). Finally the input 
presentation of (x^y) = (1,1) cancels the effect of the 
asymmetric collision of the vertical supply wave again 
deflecting the horizontal wave path {z = 0) (Fig. [6lc) 
(Mov. 3a^c [42]). 

Using this design strategy we have constructed other 
simulated gates (not, and, or, nxor) and with combi- 
nations the we have also simulated simple memory and 
arithmetic circuits (soj . 



D. Composite one bit half adder circuit 

One of the building blocks of electronic digital arith- 
metic circuits is the 1 bit half adder (ha). Two ha cir- 
cuits can be connected to make a full 1 bit adder. One 
bit adders can then be repeated connected to make an n 
bit adder. A successful construction of a HA is therefore 
demonstrative of arithmetic capabilities of this scheme. 
The HA circuit takes 2 inputs; the n bit (x) in conjunc- 
tion with a carry (y) and produces two outputs, the sum 
(S) and carry (C). A HA can be constructed from a com- 
bination of two logic gates the XOR and AND gate. There 
are two inputs {x & y) and two outputs {S & C), the 
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FIG. 6. (Colour online) Two input xOR gate (z — x^y). Left 
column sequences show wave (signal) development compared 
alongside numerical simulations on the right for the input sets 
(a) z = 0,{x,y) = [0,0]. (b) z = l,{x,y) = [0,1]|[1,0]. (c) 
z = 0, (x, y) = [1,1]- Pores of 6 SP interconnect 56 SP0 discs 
with simulation to gel projection ratio of 56 SP : 10mm 
(Mov. 3a^c 42]). 

binary sum (S) of x &l y is achieved by the XOR gate 
{S = X ®y) and inability of the configuration (overfiow) 
to present the 1 + 1 input is achieved with a carry (C) 
output, {C = X •y). 

We have conceived several simulated half adders and 
full adders circuits arranged in orthogonal [30] and hexag- 
onal \21\ networks where AND & XOR gates form building 
blocks in a larger circuit. Nevertheless it is possible to in- 
crease functional density by adding flexibility. With the 
exception of the previous diode, the above logic gates 
are constructed with discs connected with a uniform or- 
thogonal network and connection pores. More efficient 
designs are possible by permitting fiexibility in terms of 
disc size, connection pore (see diode) and inter disc con- 
nection angle. 

Figure [3 illustrates a compact HA circuit design where 
the XOR & AND gate are combined in a central single re- 
actor disc (R). Interconnecting discs (rl, r2) between the 
central reactor disc guide wave fragments to the S out- 
put to create the XOR feature (x?/(l,0) & xy{0^ 1)^5* = 
1, C = 0). The AND feature is created by the per- 
pendicular wave ejection fragments directly into the (C) 
{xy{l, 1)^S = 0, C = 1) (Mov. 4a&b [^). 

IV. DISCUSSION 

In much the same way that electronic circuits are de- 
signed and constructed we have created an alternative HA 
circuit to that presented above. In that instance the HA 
circuit is formed from the conjunction of an AND gate unit 
and a XOR gate unit. Recursive connections too have led 
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FIG. 7. (Colour online) One bit half adder circuit {S = 
C = X • y). The central reactor disc perform both 
XOR and AND logic functions to create the desired logic, (a) 
XOR component, 5* = 1, C = 0, (x, 2/) = [0, 1]|[1,0]. (b) and 
component, 5* = 0,C = 1, {x,y) = [1, 1]. Pores of 4 SP in- 
terconnect discs 56 SP0 disc, scaled by 1.2 for S and 0.8 for 
x^y with simulation to gel projection ratio of 56 SP :10mm 0. 
The circuit exploits 3 types of wave modulation, connection 
angle, disc size and pore efficacy) (Mov. 4a&b [42]). 



to the creation of memory and sequential circuits that 
can be combined with diodes and logic gates to create 
other circuits. Beyond the limitations of our own labora- 
tory resources we foresee no reason preventing successful 
migration of these additional circuits. 

Although these circuits share some similarity to their 
electronic counterparts, there are differences. Reagents, 
the chemical analog of the power supply is embedded 
and distributed throughout the circuit substrate and not 
supplied from one source. A circuit can thus activate 
for short periods in absence of any reagent refreshment. 
The vector component of the reaction wave presents an- 
other interesting contrast. Temporally separated waves 
can share a common channel without interference. It is 
therefore possible to create a 'cross-roads' junction where 
wave signals travelling in different directions can share 
the same substrate without interference ^30|. Further- 
more the optical circuit projection onto a homogeneous 
substrate permits substrate reuse with preservation cur- 
rent signals. One operating circuit projection can be 
instantaneously replace with another circuit. Waves in 
target locations could trigger circuit changes and con- 
ceivably also incorporate other existing wave signals into 
the new circuit. 

At a scale where excitation waves are visible with the 
human eye this (BZ) reaction diffusion process is rela- 
tively slow. In these experiments discs are ~ 10mm 
in diameter and waves propagate at approximately ^ 
0.5mms~^ resulting in a 205 propagation delay per disc. 
Wave front geometry is an emergent property of the 
chemical composition and substrate, that develops in- 
dependently to the initiation geometry. Reducing the 



disc physical dimensions therefore requires an propor- 
tionate increase in chemical kinetics to maintain func- 
tional equilibrium. For example; if the chemical kinetics 
could be adapted to reduce the wave front size, then a 
reduction in diameter to that of a typical neuron from 
(lOmm^O. 004mm) could at least produce a linear ad- 
justment in disc propagation delay 20/2500 = Sms; in 
comparison the refractory phase of a typical neuron is 
~ 1ms. 

This research is an exploratory component within a 
wider collaborative project that aims to create functional 
networks of lipid encapsulated BZ vesicles [26|. The lipid 
membrane and the non-linear oscillatory nature of the 
BZ medium encodes some of the features apparent in 
biological information processing; there are systematic 
analogies between electro-biochemical neurons and BZ 
discs, such as pores^ synapses and chemical—^ electrical 
excitation and refraction. 



V. FUTURE WORK 

To demonstrate the computational abilities of BZ en- 
capsulated discs we have combined Boolean information 
representation and geometric collision style manipula- 
tion. The essential components required for a univer- 
sal computer could be created with this combination. 
Nevertheless other modes of information representation 
and manipulation are possible in these networks. In an 
oscillatory mode, information could be coded in pulse 
phase relationships and manipulated by interconnected 
temporal associations, reminiscent of natural neural in- 
formation processing. Information is processed proba- 
bilistically along the edge of the BZ instability threshold 
through interconnected membrane pores. Larger reac- 
tor chambers could possibly be connected by strings of 
smaller discs or vesicles. 

Our explorations have so far been restricted in 2 dimen- 
sions. Discs have been substituted for cross sections of 
interconnected vesicles. In future studies we plan to ex- 
tend our simulations into the third dimension. We spec- 
ulate that similar modulation of wave cones could also be 
possible leading to 3D vesicle computational or adaptive 
circuits. 

Designing anything beyond all but the simplest devices 
in such schemes as is the case with neural networks, can 
be uncertain and complex. We intend to create discs ca- 
pable of adaptation in order to learn solutions as well 
as employing evolutionary strategies to search for static 
structural solutions. Currently we are exploring an evo- 
lutionary strategy to replicate the functionality of our 
manual designs of logic gates and arithmetic circuits. 
Success there opens the possibility on solving computa- 
tional tasks for which solutions are currently protracted 
in conventional systems. 
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VI. CONCLUSION 



Concomitant within our long term project aims we 
have developed a prototype information processing sys- 
tem from interconnected arrangements of BZ encapsu- 
lated discs as analogs of BZ vesicles. Expanding exci- 
tation waves sustained in a sub-excitable BZ substrate 
represent discrete quanta of information. The inter- 
connecting disc pores have a stabilising effect on the 
waves as they propagate from one disc to another. Con- 
nections, pore efficiency, connection angle and relative 
disc size are used to manipulate waves to create explicit 
chemical information processing devices. 

Networks of these units share some features apparent 
in biological information processing. Whilst inter-neuron 
communication is predominately electrical, modulation 



of that activity is chemical. In the case of individual 
neurons modulation dominates at the synaptic junction. 
The membrane pore between two vesicles can be consid- 
ered a simple analog of the synapse, a small contact area 
that can modulate signals in between vesicles. Similarly 
the electric upstate firing and downstate quiescence of 
neural signalling is an analogue of chemical excitation 
and refraction. 
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